The thermal maturity of sedimentary basins as revealed by magnetic mineralogy

The thermal evolution of sedimentary basins is usually constrained by maturity data, which is interpreted from Rock-Eval pyrolysis and vitrinite reflectance analytical results on field or boreholes samples. However, some thermal evolution models may be inaccurate due to the use of elevated maturities measured in samples collected within an undetected metamorphic contact aureole surrounding a magmatic intrusion. In this context, we investigate the maturity and magnetic mineralogy of 16 claystone samples from Disko-Svartenhuk Basin, part of the SE Baffin Bay volcanic margin. Samples were collected within thermal contact metamorphic aureoles near magma intrusions, as well as equivalent reference samples not affected by intrusions. Rock-Eval pyrolysis (T_max), and vitrinite reflectance (Ro) analysis were performed to assess the thermal maturity, which lies in the oil window when 435°C ≤ T_max ≤ 470°C and 0.6%–0.7% ≤ Ro ≤ 1.3%. In addition, we performed low- (<300K) and high-temperature (>300K) investigations of isothermal remanent magnetization to assess the magnetic mineralogy of the selected samples. The maturity results (0.37% ≤ Ro ≤ 2%, 22°C ≤ T_max ≤ 604°C) show a predominance of immature to early mature Type III organic matter, but do not reliably identify the contact aureole when compared to the reference samples. The magnetic assemblage of the immature samples consists of iron sulphide (greigite), goethite and oxidized or non-stoichiometric magnetite. The magnetic assemblage of the early mature to mature samples consists of stoichiometric magnetite and fine-grained pyrrhotite (<1 μm). These results document the disappearance of the iron sulphide (greigite) and increase in content of magnetite during normal burial. On the other hand, magnetite is interpreted to be the dominant magnetic mineral inside the contact aureole surrounding dyke/sill intrusions where palaeotemperatures indicate mature to over-mature state. Interestingly, the iron sulphide (greigite) is still detected in the contact aureole where palaeotemperatures exceeded 130°C. Therefore, the magnetic mineralogy is a sensitive method that can characterize normal burial history, as well as identify hidden metamorphic contact aureoles where the iron sulphide greigite is present at temperatures beyond its stability field.     

Structural Setting and Kinematic Analysis of the Halaban Region,Eastern Arabian Shield,Saudi Arabia

The vorticity analysis technique was applied to measure the different lithological units,such as schist,metagranite and metavolcano-sedimentary rocks,which are present in the Halaban region.This work aims to interpret the relationship between the different lithologies and the tectonic setting,in order to elucidate the nature of kinematic analysis in the Halaban region.The kinematic analyses were applied to feldspar porphyroclasts,quartz and hornblende for twentysix samples.The kinematic vorticity number (W_m) for deformed rocks in the study area ranged from～0.6 to 0.9.The direction of the long axes for finite strain data (X axes) revealed a WNW trend with shallow dipping.The direction of the short axes for finite strain data (Z axes) were represented by vertical with associated horizontal foliation.The results of the kinematic vorticity and strain analyses are characterized by simple shear with different degrees of deformation in the Halaban region.Furthermore,our finite strain data shows no significant volume change during deformation.The subhorizontal foliation was synchronized with thrusting and deformation.Furthermore,throughout the overlying nappes,the same attitudes of tectonic contacts are observable,the nappes in the orogens being formed from simple shear deformation.     

Magnetic and DC resistivity investigation for groundwater in a complex subsurface terrain

Subsurface structures associated with hard rocks are very important for groundwater. Wadi Fatima runs through the volcanic and metamorphic rocks of the Arabian Shield which are characterized by higher magnetization than the overlaying alluvium sediments. Magnetic and direct current (DC) resistivity methods have been used for groundwater exploration in the northern part of Wadi Fatima. The magnetic survey was used mainly to map the subsurface structures, using analytic signal algorithm, of the study area. The DC resistivity method was applied to describe the lithologic domain as a function of depth, depending on their electrical property contrasts where it provided a good indication for water bearing formations. The magnetic and DC resistivity interpretations were confirmed by drilling which have provided a clear idea about the hydrogeological regime of the study area. The selected drilled well is successfully productive and it produces 30 m³/h.     

On the rise of strongly tilted mantle plume tails

The rise of an initially horizontal, buoyant cylinder of fluid through a denser fluid at low Reynolds number is used to look at the ascent of strongly tilted mantle plumes through the mantle. Such ascents are characterized by (1) the growth of instabilities and (2) the development of a thermal wake downstream. Three-dimensional numerical experiments were carried out to examine these features. An hybrid particle-in-cell finite element method was used to look at the rise of non-diffusing cylinders and, a standard finite element method was used to look at the diffusing case. First the experiments show that the timescale of the fastest growing instability vary with the Rayleigh number and the viscosity ratio. In particular the growth rate decreases as the Rayleigh number decreases, in agreement with our analysis of the laboratory experiments of Kerr et al. (2008). Second the experiments show that the length of the thermal wake increases with the Rayleigh number but the change in viscosity has almost no influence on the wake length. Applied to strongly tilted mantle plumes we conclude that such plumes cannot be unstable given the plume timescales. We also discuss the application of this conclusion to weakly tilted plumes. Besides, this study allows to predict that mantle plumes are unlikely to have developed a significant thermal wake by the time they reach the surface. Finally, the resolution that is required to allow for the growth of mantle plume tails by combined diffusion and thermal entrainment is shown to represent a challenge for the large scale mantle convection simulations.     

Petrogenetic modeling of rock variety in the Khalkhab–Neshveh pluton,NW of Saveh,Iran

The Khalkhab–Neshveh (KN) pluton is a part of Urumieh–Dokhtar Magmatic Arc and was intruded into a covering of basalt and andesite of Eocene to early Miocene age. It is a medium to high‐K, metaluminous and I‐type pluton ranging in composition from quartz monzogabbro, through quartz monzodiorite, granodiorite, and granite. The KN rocks show subtle differentiation trends strongly controlled by clinopyroxene, plagioclase, hornblende, apatite, and titanite, where most major elements (except K₂O) are negatively correlated with SiO₂; and Al₂O₃, Na₂O, Sr, Eu, and Y define curvilinear trends. Considering three processes of magmatic differentiation including mixing and/or mingling between basaltic and dacitic magmas, gravitational fractional crystallization and in situ crystallization revealed that the latter is the most likely process for the evolution of KN magma. This is supported by the occurrence of all rock types at the same level, the lack of mafic enclaves in the granitoid rocks, the curvilinear trends of Na₂O, Sr, and Eu, and the constant ratios of (⁸⁷Sr/⁸⁶Sr)_i from quartz monzodiorite to granite (0.70475 and 0.70471, respectively). In situ crystallization took place via accumulation of plagioclase and clinopyroxene phenocrysts and concentration of these phases in the quartz monzogabbro and quartz monzodiorite at the margins of the intrusion at T ≥ 1050°C, and by filter pressing and fractionation of hornblende, plagioclase, and later biotite in the granitoids at T = ～880°C.     

Impact of model complexity and multi-scale data integration on the estimation of hydrogeological parameters in a dual-porosity aquifer

This study investigates the impact of model complexity and multi-scale prior hydrogeological data on the interpretation of pumping test data in a dual-porosity aquifer (the Chalk aquifer in England, UK). In order to characterize the hydrogeological properties, different approaches ranging from a traditional analytical solution (Theis approach) to more sophisticated numerical models with automatically calibrated input parameters are applied. Comparisons of results from the different approaches show that neither traditional analytical solutions nor a numerical model assuming a homogenous and isotropic aquifer can adequately explain the observed drawdowns. A better reproduction of the observed drawdowns in all seven monitoring locations is instead achieved when medium and local-scale prior information about the vertical hydraulic conductivity (K) distribution is used to constrain the model calibration process. In particular, the integration of medium-scale vertical K variations based on flowmeter measurements lead to an improvement in the goodness-of-fit of the simulated drawdowns of about 30%. Further improvements (up to 70%) were observed when a simple upscaling approach was used to integrate small-scale K data to constrain the automatic calibration process of the numerical model. Although the analysis focuses on a specific case study, these results provide insights about the representativeness of the estimates of hydrogeological properties based on different interpretations of pumping test data, and promote the integration of multi-scale data for the characterization of heterogeneous aquifers in complex hydrogeological settings.     

PROMETHEE II and fuzzy AHP: an enhanced GIS-based landslide susceptibility mapping

This paper describes the application of a well-known multi-criteria decision-making technique, called preference ranking organization method for enrichment evaluation (PROMETHEE II), in combination with fuzzy analytical hierarchy process (FAHP), as a weighting technique to explore landslide susceptibility mapping (LSM). To this end, eight landslide-related geodata layers of the Minoo Dasht located in the Gorgan province of Iran, involving slope, aspect, distance to river, drainage density, distance to fault, mean annual rainfall, distance to road and lithology have been integrated using the PROMETHEE II enhanced by FAHP technique. Afterward, the receiver operating characteristics (ROC) curves for the proposed LSM were drawn using an inventory of landslides containing 83 recent and historic landslide points, and the area under curve = 0.752 value was calculated accordingly. Additionally, to further verify the practicality of such susceptibility map, it was also evaluated against the landslide inventory using simple overlay. The outcome was that about 11 % of the occurred landslide points fall into the very high susceptibility class of the LSM, but approximately 52 % of them indeed fall into the high and very high susceptibility zones together. Also, it resulted that no recorded landslide occurred in the zone of very low susceptibility. According to the results of the ROC curves analysis and simple overlay evaluation, the produced map has exhibited good performance.     

Evaluation of distributed recharge in an upland semi-arid karst system: the West Bank Mountain Aquifer,Middle East

Assessment of recharge in a structurally complex upland karst limestone aquifer situated in a semi-arid environment is difficult. Resort to surrogate indicators such as measurement of spring outflow and borehole discharge, is a common alternative, and attempts to apply conventional soil moisture deficit analysis may not adequately account for the intermittent spate conditions that arise in such environments. A modelling approach has been made using the West Bank Mountain Aquifer system in the Middle East as a trial. The model uses object oriented software which allows various objects to be switched on and off. Each of the main recharge processes identified in the West Bank is incorporated. The model allows either conventional soil moisture deficit analysis calculations or wetting threshold calculations to be made as appropriate, and accommodates both direct recharge and secondary recharge. Daily time steps enable recharge and runoff routing to be calculated for each node. Model runs have enabled a series of simulations for each of the three aquifer basins in the West Bank and for the whole of the West Bank. These provide recharge estimates comparable to those prepared by earlier workers by conventional means. The model is adaptable and has been successfully used in other environments.     

Artificial neural network inversion of magnetic anomalies caused by 2D fault structures

A new approach is proposed to interpret magnetic anomalies caused by 2D fault structures. This approach is based on the artificial neural network inversion, utilizing particularly modular neural network algorithm. The inversion process is implemented to estimate the parameters of 2D fault structures where it has been verified first on synthetic models. The results of the inversion show that the parameters derived from the inversion agree well with the true ones. The analysis of noise has been studied in order to investigate the stability of the approach where it has been tested for contaminated anomalies with 5 and 10 % of white Gaussian noise. The results of the inversion provide satisfactory results even with contaminated signals.The validity of the approach has been demonstrated through real data taken from New South Wales, Australia. A comparable and satisfactory agreement is shown between the inversion results of the neural network and those from techniques published in literatures.     

Radiative and free-convection effects on the oscillatory flow past a vertical plate

The interaction of free convection with thermal radiation of the oscillatory flow past a vertical plate is studied. The Rosseland approximation is used to describe the radiative heat flux in the energy equation.